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We have recalculated the Mikaelian and Smith radiative corrections to the Dalitz decay 7r° —>■ 
e''"e “7 beyond the soft-photon approximation, i.e. over the whole range of the Dalitz plot and with 
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one-loop level and the virtual muon loop contribution. The results can then be used also for heavier 
particles in the final state. 
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I. INTRODUCTION 

Right after the process 7r° —>• 77, the second most 
important decay channel of a neutral pion is the 
Dalitz decay tt® —>■ e“*"e“7 with a branching ratio 
(1.174 ± 0.035)% [I|. This decay was named after 
Richard H. Dalitz, who first studied it in Ref. 0. Experi¬ 
mental data of this process provide information about the 
semi-off-shell pion transition form factor J7«77* (Q^/M^) 
in the timelike region and in particular its slope param¬ 
eter a. 

Radiative corrections to the total decay rate of the 
Dalitz decw 7r° —>■ e+e“7 were first addressed by 
D. Joseph [^. The pioneering study of the corrections 
to the differential decay rate was done by B. E. Lautrup 
and J. Smith using the soft-photon approximation. 
This analysis was soon after extended by K. O. Mikaelian 
and J. Smith by hard-photon corrections to the whole 
range of the bremsstrahlung photon energy. As one of the 
main results of their work the table of radiative correc¬ 
tions S{x,y) to the leading-order (LO) differential decay 
rate was presented. 

It turned out that such a table would be very useful for 
the Monte Carlo simulations in experiments covering 7r° 
decays, e.g. the NA48 experiment at CERN Q. In prac¬ 
tice, for the table of values 5{x,y), which was published 
in Ref. [^, an interpolation or extrapolation procedure 
needs to be used in order to get the radiative correction 
at any desired point of the Dalitz plot. This might lead 
to a large uncertainty. 

We have therefore recalculated, generalized and ex¬ 
tended the results presented in Ref. [Sj and prepared the 
code which can give a value at any kinematical point 
{x,y). As we have not neglected the higher-order terms 
in the electron mass and included also the muon loop con¬ 
tribution to the vacuum polarization insertion correction, 
our result can be in principle also applied to the other 
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related processes. The decay of an eta meson to a muon 
pair and a photon, where the masses of the final-state 
particles are not anymore negligible in comparison to the 
decaying pseudoscalar, is such an example. On the other 
hand, when an eta meson and its decays come into play, 
some peculiarities inevitably appear. We comment on 
this a little in the present work but postpone the details 
and the results of the radiative corrections for this case 
to the paper in preparation. Nevertheless, we try to be 
as general as possible considering the presented results so 
one can utilize the formulas without modifications later 
on. 

To proceed even further we have also included the one- 
loop one-photon-irreducible contribution, which was con¬ 
sidered to be negligible in the original paper @ due to 
its proportionality to the Iraton mass. This statement 
had been corrected in Ref. [3 many years before the de¬ 
bate about this issue was closed; see e.g. Refs. @,0. We 
provide here a complete calculation of this contribution 
making no approximations considering the lepton masses 
and energy of the photon. We show that this correction 
is indeed important and changes significantly the values 
of entries stated in Table I of Ref. @ especially for a large 
invariant dilepton mass. 

Let us also mention that a systematic treatment of the 
next-to-leading-order (NLO) corrections to the Dalitz de¬ 
cay of a neutral pion in the framework of chiral pertur¬ 
bation theory with dynamical leptons and photons was 
studied in Ref. M- Here we will also use some results of 
this work. 

It is worth it to notice that throughout the paper we 
stick to the notation which was used in Ref. @ using only 
minor modifications. Even though some of the names 
may appear to be clumsy, we believe that it would be 
confusing to do otherwise. Naturally, such an approach 
is also very convenient for the reader who is familiar with 
the original work. 

Our paper is organized as follows. We recapitulate 
first some basic facts about the LO differential decay 
width calculation in Sec. im Then we proceed to the re¬ 
view of the NLO radiative corrections in the QED sec¬ 
tor in Secs. MMl and El In particular, in Sec. imi we 
discuss the virtual corrections including the muon loop 
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Figure 1. Leading order diagram of the Dalitz decay 7r° ^ 
e^e ~7 in the QED expansion. 

contribution, in Sec. IIVI we introduce the one-photon- 
irreducible contribution and in Sec. El we describe the 
bremsstrahlung correction calculation. Some technical 
details together with extensive results concerning the 
bremsstrahlung contribution to the NLO correction have 
been moved to the Appendixes. 


II. LEADING ORDER 


First, let us briefly introduce some basic notation. In 
what follows we denote the four-momenta of the neutral 
pion (of the mass M), electron (mass to), positron and 
photon by P, p, q and fc, respectively. We also introduce 
common kinematic variables x and y defined as 


{P + ^_2 P ■ ip-q) 

M2 ’ ^ M2 (1 - x) 


( 1 ) 


where a: is a normalized square of the total energy of the 
e“'"e“ pair in their center-of-mass system (CMS), or sim¬ 
ply of the electron-positron pair invariant mass. The vari¬ 
able y has then the meaning of the rescaled cosine of the 
angle between the directions of the outgoing photon and 
positron in the e+e" CMS. If we introduce v = 2m/M 
and 


II 

II 


(2) 

we can write the limits on x 

and y as 


X e [ 1 ^ 2 , 1 ] , 

y e [-/3,/3]. 

(3) 


The leading-order diagram of the Dalitz decay 7r° —>• 
e“'"e “7 is shown in Fig. [TJ The shaded blob corresponds 
to the neutral pion semi-off-shell transition form factoi0 




^ As it also follows from the definition 0. we will use shortly 
J^(0) = (0), which is complementary to the dou¬ 

bly off-shell transition form factor taken at the photon point 
(0) = *7^7^07*7* (0, 0) = which at the LO of the 

chiral expansion is equal to = — l/(47r^i^). 


which is related to the doubly off-shell tran- 
sition form factor JVo 7 * 7 * (Qi/AI^, = 

defined as 


J (5) 

= (/2/m2, (P - 0 VM^) , 

by P,o^^^.(Q2/m 2) = P,o^.^.(0,QVM2). In Eq. 0, 
/ is a dimensionless function, which can be linearly ex¬ 
panded in the chiral perturbation theory in terms of the 
slope parameter a as follows 


/(2;)~I-haz. (6) 

In our case it then holds = {P — k)^ = M'^x and for 
the leading-order matrix element in the QED expansion 
we can write 

^M^°{p,q,k) = ^P{x)e*Pik) 

X 2to [u{p, m)jp^j5v{q, to)] 

^ (7) 

+ [Hp, w) [7p (k-p)- ppli:] 75u(g, TO)j 

- u{p, to) [ 7 p {k- q)- gp|f] -i5v{q, TO)j | . 


Summing the modulus squared of the previous result over 
the fermion spins and photon polarizations and taking 
into account that, in general, in terms of variables x and 
y it holds 


dr(a;, 2 /) = —^|7W(a;,?/)|^(I-a;)da;dy, ( 8 ) 


the differential decay rate then reads 


da:d?/ ( 871 )^ 2 


(l-xf 




l + y^ + - 
(l-x) 


Ax 


1 + r 



(9) 


Here we have used the LO expression for the decay rate 
of the neutral pion main decay mode 

T n , / m9 , n 

( 10 ) 

Integrating m over y we find 


dF^O 

dx 


©i/(.)P 


^LO 

77^^77 


8/3 (1 - x) 


3 r 


4x 


, 2 1 


1 -b 


2 x 


( 11 ) 


Moving beyond the leading order, it is convenient to 
introduce the NLO correction <5 to the LO differential 
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decay width, which can be in general defined as (in the For the correction stemming from the vacuum polariza- 
case of the two-fold differential decay width) tion insertion in Fig. HJa) we can write 


5{x,y) 


^2pNLO /^2pLO 

dx&y / AxAy 


or (in the one-fold differential case) 


(5(a:) 


dpNLO /dr^o 
da; / dx 


( 12 ) 


(13) 


Such a correction can be divided into three parts empha¬ 
sizing its origin 

5 -p -P (5®® . (14) 

Here, stands for the virtual radiative corrections, 
^i 7 iR Pqj, .j-dg one-photon-irreducible contribution, which 
is treated separately from in our approach, and (5®® 
for the bremsstrahlung. Having knowledge of d{x,y), we 
can calculate S{x) as a trivial consequence of previous 
equations using the prescription 


5(x) 


3 1 

^(1 + i) 


S{x,y) 


, 2 -| 


1 + y + 


dy. (15) 


In the following sections, we discuss the individual con¬ 
tributions one by one. 


III. VIRTUAL RADIATIVE CORRECTIONS 




(a) (b) 

Figure 2. Virtual radiative corrections for 7r° —>■ e'''e “7 
process; vacuum polarization insertion (a) and correction to 
the QED vertex (b). 


From the interference terms of the LO diagram shown 
in Fig. [1] with the one-loop diagrams presented in Fig. [51 
we get NLO virtual radiative corrections which can be 
written as @ 

6'''’^^{x,y) = 2Re n(x) -P Fi{x) + 

(16) 

or (through the formula (1151) 1 as 
<5™‘(x)=2Re|-n(x)+Fi(x) + ^^^^| . (17) 


2F2{x) \ 

l + y^ + ^J 


n(x) = ne(x)-pn^(x). (is) 


Here we have explicitly written not only the contribution 
coming from the electron loop as it was done in Ref. 
but also from the muon loop. This becomes both nec¬ 
essary and convenient when one goes beyond the decay 
7r° —>■ e+e^y, which we discuss throughout this text, and 
proceeds to the process rj —>■ pt+/r“ 7 . In the end, we 
might then simply do the exchange rrie -O- (we and 
stand for the electron and muon mass, respectively) 
in the expression for the correction 5 to vary the final- 
state lepton masses. Let us remark that independently 
of the considered processes, the loop with the lightest 
fermion is of the greatest importance. Thus, taking only 
the electron loop into account (i.e. leaving the muon part 
in (1181) 1 and performing simply the tempting lepton mass 
substitution in the whole expression, we would miss out a 
very important contribution. Obviously, the vacuum po¬ 
larization insertion defined in a way shown in (I18p stays 
after such an operation intact, as desired. The other op¬ 
tion would be to treat separately the final-state lepton 
masses m and the masses of the particles in the vacuum 
polarization insertion loops me and m^. This more uni¬ 
versal approach was used in the code which comes with 
the paper. Let us now introduce for the later convenience 


7 = 7 ( 2 ;) 


1 - I3{x) 

1 -I- j3{x) ' 


(19) 


The individual terms used in (1181) are then defined as 


a 

1 



f 

TT 

~9 

3 

^ 2 xy 

V 


( 20 ) 

In the above formula, I stands for e or /r in the loop and 
changes the meaning of the so far used electron mass m 
in the definitions of v, j3 and 7 to me or m^. Unlike 
in Ref. @ where only the real part of (l20l) above the 
threshold x = vf is shown, we quote here the full expres¬ 
sion valid in all kinematical regimes. This is necessary 
to get right the contribution from the charged fermion 
loop when the transferred momentum is not sufficiently 
large to produce the real pair, i.e. for x < vj, and lacks 
therefore the imaginary part. This situation for instance 
appears (at least for a part of the kinematical region) 
when the pseudoscalar decays to the electron-positron 
pair via the muon loop. For the purpose of real algebra 
used in the code (i.e. to avoid complex logarithms and so 
on) we can extract the real part of (|20)) . For an arbitrary 
mass of the charged loop fermion we find 


Re{/ 3 flog[- 7 ^]} 

= —2|/3^| ■|0(/3|) arctanh/3^ -|- 6{—Pf) arctan|^| . 

( 21 ) 
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In the following, we stick exclusively back to the process 
7 r° —>■ and m then denotes the outgoing electron 

mass as before. Finally, for the electromagnetic form 
factors Fi{x) and F 2 {x) stemming from the QED vertex 
correction in Fig. [2jb) we have 


, , a ( 1 + 2/32 


1 + P^ 
2/3 


log(- 7 ) 


1 


Li 2 (l- 7 ) + -log (- 7 ) 




1 + 


I +/32 

2/3 


log(- 7 ) 


m 

logy 


( 22 ) 


and 


F2{x) 


-^log(- 7 ). 
TT Axp 


(23) 


In the above formulas, Li 2 stands for the dilogarithm and 
A is the infrared cutoff. To extract the real parts from the 
previous terms (1221) and (l23l) (in a sense of applying the 
operator Re), in the kinematically allowed region where 
M^x > 4to 2 we use log(— 7 ) = log( 7 ) + itt, since 0 < 
7 < 1. Thus it is straightforward to see that the real 
parts of Fi(x) and ^ 2 ):^ indeed coincide with the form 
factors stated in Ref. Q including the Coulomb term 
proportional to 


IV. ONE-PHOTON-IRREDUCIBLE VIRTUAL 
RADIATIVE CORRECTION 


+ cross + CT 

(a) 



Figure 3. One-loop one-photon-irreducible contribution 
LO Feynman diagrams for 7r° e'*'e “7 process considering 

the QED and yPT expansion: triangle diagrams and related 
counterterms (a) and box diagram (b). Note that “cross” 
accounts for a diagram with a photon coming from the out¬ 
going positron line. “CT” then stands for two counterterm 
diagrams necessary to compensate the UV divergent parts of 
the related Feynman diagrams. 



Here we will summarize the most important results which 
are necessary to proceed toward our purpose considering 
NLO corrections to the decay 7 r° —?► e+e^y. Let us em¬ 
phasize that this contribution was not included in the 
calculations performed in Ref. On the other hand, it 
was shown later on in Ref. Q within the limit to —>■ 0 
that there is no point in treating the I 7 IR correction as 
negligible. In the following we show the results of the 
calculation beyond this massless limit. For the reasons 
specified in the previous sentences we have devoted to 
this contribution a separate section, even though it is of 
course just one additional virtual radiative correction. 

Until now we have not considered any particular form 
of the semi-off-shell form factor F{x) in our calculations. 
To get the one-photon-irreducible contribution in a closed 
form, it is though necessary to choose a concrete form of 
Accordingly, we should consider at this mo¬ 
ment a general doubly off-shell pion transition form fac¬ 
tor (Z 2 /m 2 , (P — /) 2 /m 2 ], where I stands for a 

loop momentum. In Fig. |3 we can see the LO of the 
considered contribution in chiral perturbation theory. In 
such a limit we take the constant = — I/( 47 r 2 F) 

as the local form factor and it is thus clear due to the 
power counting that a counterterm is needed. The fi¬ 
nite part of such a counterterm renormalized at scale /i 
is governed by the parameter x*-‘^^(/i), which corresponds 
to the high-energetic behavior of the complete form fac¬ 
tor. This can be theoretically modeled e.g. by the lowest- 
meson-dominance (LMD) approximation to the large-Vc 
spectrum of vector-meson resonances yielding the value 
= 2.2 ± 0.9 [ 13 , which can be further used 
for numerical results. The dependance of the correction 
15 I 7 IR Qjj ^(r) neglected for the values given by rel¬ 

evant models as well as experiments, when a decay with 
the electrons in the final state is taken into account. We 
will comment on this in the end of this section. Let us em¬ 
phasize in a more straightforward way that using “only” 
the LO expansion of the form factor is compensated by 
the effective value which differs for particular mod¬ 
els. One gets the model corresponding value of for 
instance from the matching to the full calculation. In 
this sense one loses no information. On the contrary, the 
model dependence of any such result can be conveniently 
altered easily just by changing the value of ■ 

The total matrix element covering all the diagrams rep¬ 
resented in Fig. 13 can be written in such a form which 
manifestly satisfies the Ward identities for the conserved 
electromagnetic vector current 


iMi-tmip, g, k) = - 


X [{k-p)qp - {k ■ q)Pp] [u(p,TO)75i;(g,TO)] 

+ A(x, y) u{p, to) [7p {k-p)- ppl/;] Xbv[q, m) 


( 24 ) 


One-photon-irreducible (I 7 IR) contributions were ex¬ 
tensively studied in Ref. [llj in connection with the 
bremsstrahlung correction to the 7 r° —>■ e“'"e“ process. 


- A{x, -y) u{p, to) [ 7 p ik-q)- qp^] 75 ^( 9 , to) 
+ T{x, y) [u{p, m)xpi(:x 5 v{q, to)] | . 
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Here P, A and T are scalar form factors, the explicit form 
of which can be found in Appendix A of Ref. [TJ . 

To get the NLO one-photon-irreducible part of the cor¬ 
rection 6 we need to consider the interference term of LO 
matrix element © and the I 7 IR contribution (l24l) and 
sum it over the photon polarizations with the result 


■^ 17 Ir( 2 ;: 2 /) = fc)]*-Adi7iR(p,g, fc) 


polar. 

-^- - - ^ {x)J^^o^^^4:i^T{x,y) 

A{x, y)M[x{l - yf - v'^] + (y ^ -y) |. 


(25) 


Putting the above formula into ([ 8 ]) and (IT^ and normal¬ 
izing to the LO two-fold differential decay width m we 
get finally 




dccdy 


= 2 Re <^ - 


aJ'^O(O) in^M 


-b 


A{x,y)M[x{l - yf - v'^] + {y ^ -y)j ||. 


(26) 


For our purpose we can safely set P{x) ~ J'^°(0) in the 
previous formula, considering only the leading order of 
the chiral expansion; see also m assuming the slope a 
is small. It should be mentioned, though, that such an 
approximation is only reasonable for the Dalitz decay of a 
neutral pion. For the decays of an eta meson, one should 
be more cautious and use a better treatment of the full 
form factor. 

Similarly, the dependence on the parameter cannot 
be neglected when i/ becomes significant. Indeed, con¬ 
sidering the full expression (A.5) from for the form 
factor T{x,y), one gets for the y-dependent contribution 
to (5^'>'®(a;, j/) from (l26l) 



aJ'i^O(O) 4:/2^W(/x) 1 

TT P{x) (1 - x){l - 2 / 2 ) [1 + y 2 • 


(27) 

Thus, e.g. for the decay ij —>■ /r+/i 7 the one-photon- 
irreducible contribution may be considerably model- 
dependent. This is, however, not the case for the pro¬ 
cess 7 r° —>■ e^e “7 where the contribution given in (1271) is 
suppressed in comparison to the other terms in (I26p . 


V. BREMSSTRAHLUNG 

In this section we recapitulate the approach used in 
Ref. for the bremsstrahlung correction calculation. 
We think it is useful and convenient to rewrite the whole 
story in a more detailed way so it is transparent and easily 



+ cross 

Figure 4. Bremsstrahlung corrections for 7r° —>■ e'’'e “7 pro¬ 
cess. Needless to say, “cross” stands for the diagrams with 
outgoing photons interchanged. 


understood. As usual, one can then build on that when 
a few more pieces come into play. In the Appendixes we 
then provide the results themselves. Note also that espe¬ 
cially in this section we restrict ourselves to the original 
notation used in the work Q. 

The diagrams which contribute to the Dalitz decay 
bremsstrahlung and are thus important to cancel the IR 
divergences stemming from the virtual corrections dis¬ 
cussed in Sec, mil are shown in Fig. 3 ] The corresponding 
invariant matrix element (including cross terms) can be 
written in the form 


zMbs = u{p) [r"(fc, 1) + k)\v{q)e*p{k)C{l ), (28) 

whercQ 




2 ^{i+p+q){k)m 

{l+pPqY 


' (/+^ + m) 

21 ■ p Pie 


. M + j-nn) 


21-q 


le 


(29) 


The form factor F{{1 P p P qY/M'^) can be expanded 
(assuming a is small) in the following way 




1 H" d 


(p + q) 

M 2 


(30) 


Thus for the process 7 r° —> e^e~^ it can be approximated 
by P{x), taking into account only the leading order in the 
chiral expansion. Let us also introduce Tr for the rescaled 
matrix element squared and summed over all spins and 
polarizations of final states by the relation 


|AIbsP= ^ |AlBSp = ^|.F(a;)pTr. (31) 

sp., pol. 

Inasmuch as an additional photon comes into play it is 
convenient to introduce a new kinematic variable which 


^ We use the shorthand notation for the product of the Levi-Civita 
tensor and four-momenta in which . 
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describes the normalized invariant mass squared of the 
two photons 


(fc + 0 " 

~ M 2 


(32) 


It has the similar meaning as x in the case of the electron- 
positron pair. The limits on x^ are 


M 2 


< Xj < x 


max 

7 


1 + X — \ ix + ^{1 — xy . (33) 


The contribution of the bremsstrahlung to the next- 
to-leading order can be described (according to (fT^ l by 
the correction 


If we follow the notation of Ref. @ , we define the prop¬ 
agator denominators in the following way (suppressing 
+ie part for now) 

A = l-q, B = l-p, C = k-q, D = k-p, 

E = {p + q + l)^ , F = {p + q + ky . 

Not only is the whole amplitude invariant under the in¬ 
terchange of the two photons (and thus of k and I in 
(l28l) l. but also the operator J possesses the same sym¬ 
metry which can be written for an arbitrary function of 
the propagator denominators m as 

J[fiA, B, C, D, E, F)] = J[f{C, D, A, B, F, E)] . (42) 


in which in 

d^r^s _ 

dxdy 




d2rBS /d2rLO 


da:d?/ / dxdy 
agreement with Q we can write 
(1-x) tt^M^ 


4M(27r)8 16 

|/(x)p /a\2 

Vtt 


J 




64 


T: 


LO 




+77 


|A4bs 
( 1 - x) I J[Tr] 


(34) 


dx-, 


(35) 


The above used operator J is defined for an arbitrary 
invariant /(fc, 1) of the momenta k and I as follows 

1 r u3l j3/ 

= ^ J j^^f{k,l)S^^\P-p-q-k-l). 

(36) 

Finally, putting the LO differential decay width expres¬ 
sion ® and the previous result (1551) into (IMl) we get 


1 /a\ 4x J J[Tr] dx.y 
64 V TT/ (1 — x)2 1^1 -(- y2 _j_ 


(37) 


In the remaining part of this section we discuss the way 
the integral / J[Tr] dx-y is treated. Most of the explicit 
formulas are then moved to the Appendixes. 

Being on shell (fc2 = 0 = and in the diphoton 
center-of-mass system where P —p — q = 0(= k + l = r), 
we find 


J[f{k,i)] ^ ^ I dn f{k,k) = ^j dnjiU). 

(38) 


Here, we have used I to mark the four-momentum I 
with the opposite momentum direction, i.e. whenever 
I = {lo, 1), then I = {Iq, —1). We can come back to the 
invariant form in a known way through 

2lo = ko + lo ^{k + iy = (39) 

or, for example, due to 


The interchange of p and q (which is also a relevant sym¬ 
metry in our case) must be compensated on the level of 
the operator J by changing the y sign, thus 


J[/(A, H, C, D, E, E)] = J[f{B, A, D, C, A, F)] 



There are also some useful identities which follow from 
the definitions (HD) such as 


E-2A-2B = F-2C-2D = M^x , (44) 

F + 2A + 2B = E + 2C + 2D = M^{1-Xj) (45) 

and 

A + C = ^[{1- x){l+y)-Xj] , (46) 

M^ 

B + D = — [{l-x)il-y)-Xj] , (47) 

= iA + C) 

y^-y 

E + F = M^{l+x-xj). (48) 


It is convenient to know the above relations for two rea¬ 
sons. First, we see that we can simply trade one of the 
above defined variables for the others and thus only two 
more independent variables in addition to x, y and x^ 
(e.g. A and B) are necessary to describe the kinematics 
of our decay. On the other hand, we realize that some 
special combinations of the variables A, ..., F are invari¬ 
ant with respect to the acting of the operator J (i.e. they 
depend only on x, y and x-y). We can also combine the 
previous formulas to get some other J-invariant combi¬ 
nations. If we consider, for example, that 


A={A + C)-^[{E + 2C + 2D)-E-2D], (49) 

we find 
F' 

— - A +D =—[l + x - x^ - y{l-x)]. (50) 


_ (fcp -I- lo)po (f=o) {k + l) -p 
(ko + lo) ~ 


( 40 ) 


Such expressions are useful when we want to reduce 
the complicated J terms, arising naturally during the 
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calculation of the invariant matrix element squared, to 
the basic ones which are simple to handle. First, we 
use the above stated relations to simplify the numerators 
(e.g. we get rid of ^ in a term like A/{DE) using the 
relation (ISOt lFI Then also the denominators are treated. 
For example, consider the term J[1/{ACEF)]. Then 

1 1 1 {A + C){E + F) 

AC EE ~ {A + C){E + F) ACEF 

1 / 1 1 1 1 \ 

“ {A + C){E + F) \AE ^ CE CjE) ' 

After applying the operator J and using the symmetry 
(H^ . we find 


J 


1 


ACEF 


{A + C)iE + F) 


J 


AE 


+ J 


CE 


(52) 


All necessary reductions of this type are summarized in 
Appendix m except for such terms which one can get 
using the discussed symmetries and (H51) . The com¬ 
putational methods used to calculate the basic terms are 
introduced in Appendix [C] For the list of the results for 
these integrals see Appendix |DJ Flere, in comparison to 
Ref. we include also the new term J[1/{A‘^E^)] which 
appears due to the fact that 0{iy'^) terms were not ne¬ 
glected in our approach. The completely reduced rescaled 
matrix element squared Tr, which represents in terms of 
J[Tr] an important ingredient for the bremsstrahlung cor¬ 
rection S^^{x,y) (cf. (l37l) '). is presented in Appendix lAl 
We believe we provide here the results in a more refined 
way in comparison with Ref. Q. 

The last step is the integration over Xj. There are 
basic integrals which behave like 1 /x-y and are divergent 
when this integration is performed if no Xy appears in 
the numerator to compensate it. The essential divergent 
integrals are J[l/A^] and J[1/{AB)]. The divergent part 
of integrals like J[l/{A^E)] and J[1/(A^£’^)] can then 
be written in terms of these essential ones. For example, 
using dm) we get 

1 _ 1 / 1 2 2B 

“ M2x ~~AE~ 



Needless to say, there are also A ^ B counterparts of 
the mentioned integrals. This unwelcome behavior can 
be extracted from the Tr expression to get the conver¬ 
gent part Trc, which can be treated numerically, and the 
divergent part Ti'd, which should be treated analytically. 
In the former case we can set A —>■ 0 and the lower bound 


^ The combination llsoll is of course in some minimalistic sense re¬ 
dundant for the considered procedure, since we can always make 
two-step substitution instead. In such a case, we would trade A 
for C using Hell and then C for E and D using gsj. 


on Xy is then zero. In the latter case the cutoff A has to 
be preserved. 

Finally, as expected, the sum of the divergent part of 
the bremsstrahlung correction y), the explicit form 

of which can be found in dA4|) , and the divergent part of 
virtual correction d™^{x,y), represented in the follow¬ 
ing formula by the electromagnetic form factor F'i(x), in 
particular 


+ 2Re{Fi(x)} , (54) 

is IR finite. In other words, terms proportional to 
log to/A cancel each other in the final formula of the cor¬ 
rection S(x, y). 

In the end of this section, let us go back to Eq. (1501) . In 
cases when the slope a is no longer negligible in compar¬ 
ison to I, one should consider the entire right-hand side 
of (l30l) instead of only T{x) alone. It is then necessary 
to go beyond the approach used in Ref. @ . If we square 
the bremsstrahlung matrix element (I28L we get for the 
simple case with a = 0 which we have treated so far 

I-MbIT = \h,i + Ak? = \h,i?E\h,k?+2IliIi,k • (55) 
Here we have denoted 


Ik,i = l)v{q)e*pik)el{l) (56) 

and likewise for Using the building blocks of the “no¬ 
slope” matrix element modulus squared (l55l) and consid¬ 
ering the expansion (1301) we find the correction for the 
bremsstrahlung expression 


|M“f|2 = [l + a(I-x-x,)]|M“iO|2 

+ {JEW-W) . 


(57) 


If we apply the operator J and take into account the 
symmetry (|42l) . the previous formula can be boiled down 
to 


J 


iMitr - \M 


BS 


= 2a< (I — X — Xy) J 


T* T 


AJ 


{A + B) 


M2 


\Ik.l 


(58) 


This expression can be calculated along the same lines 
as T[Tr]. One then gets a similar expression to Trc in 
(|A2p including some new integrals. These need to be 
calculated in addition to the known basic terms. Note 
that there is no divergent part in (1551) which needs to be 
treated separately. 

The above correction does not need to be considered in 
the decay 7 r° —>■ e+e “7 so we do not present the related 
results in this paper. On the other hand, it becomes 
important when treating the eta meson decays. 



































VI. RESULTS 


For the reader’s convenience, we put here together the 
individual pieces (ITBl) , (1^ and (1571) and write the overall 
NLO correction 


S(x, y) = y) + y) + y) 

2F2[x) 


= 2 Re ~n(x) + Fi {x) + 
M 




A(x,y)M[T(l - yf - + {y ^-y) || 


1 /a\ Ax f J[Trc]da;.y 


64 


© 


TtJ (1-xy [l+y2_^id] 


+ Sjj^(x,y). 


(59) 


Here, the convergent part of the rescaled bremsstrahlung 
invariant matrix element squared (to be integrated over 
Xj numerically) Tr^ is given by (IA2I) and the analytically 
integrated divergent part of the bremsstrahlung correc¬ 
tion ( 5 q®(x, 2 /) is shown in (IA4I) . Let us recall that the 
explicit formulas for the scalar form factors A and T can 
be found in Appendix A of Ref. [HI- 

Taking the result (|59)l and using the formula m, 
we get the overall correction to the one-fold differential 
leading-order decay width, which is shown in Fig. [5j For 
comparison, also the sum 5™*(a:) -I- (5®®(a:), which would 
have corresponded to the correction presented in the orig¬ 
inal paper if the terms and the muon loop had 

not been omitted, and one-photon-irreducible contribu¬ 
tion are shown. We see that in the case of the 

decay 7 r° —>■ the I 7 IR correction is negative for 

the whole range of values of x and enhances thus the ef¬ 
fect of the sum (a:) -I- <5®® (a;) which is also negative in 
a wide range of x. 

Taking into account all the discussed contributions, a 
similar table of values of correction 5{x, y) as it was pro¬ 
vided in the original work Q, can be produced at the 
very same points according to (1551) ; see Table HI Consid¬ 
ering the contributions introduced in this work but left 
out in Ref. Q, the I 7 IR correction is the most impor¬ 
tant one, especially for large x. The correction of the 
old Mikaelian and Smith values is significant and greater 
than 10% already for x ~ 0.5. This can be visible in 
Fig. [5] and also from the difference of the entry values 
between the Table I in Ref. @ and Table H] in the present 
work, provided the remaining contributions are not sig¬ 
nificant. Indeed, the muon loop vacuum polarization in¬ 
sertion contribution, which is independent on y, grows 
nearly linearly with x from i5)))ioop(0.01, y) = —0.0005% 
up to <^)))ioop(0.99, y) = —0.0616% and is thus negligible. 
A similar conclusion holds then also for the 0{v^) con¬ 
tribution, which is most significant for small x with the 
value 5(^p(0.01,0) = 0.0035%. 

With our present knowledge, we are now in a position 
to calculate the correction to the integrated decay width. 



Figure 5. The overall NLO correction 5{x) for the decay tt® —>■ 
e^e “7 calculated according to the formula (1591) (solid line) in 
comparison to its constituents. The sum -I- <5®® (a:) is 

depicted as a dashed line and the one-photon-irreducible con¬ 
tribution is shown as a dash-dot line. The divergent 

behavior of S{x) near a; = ~ 0 has the origin in the electro¬ 

magnetic form factor Fi {x) and is connected to the Coulomb 
self-interaction of the dilepton at the threshold. 


In this case, the transition form factor F{x) cannot be 
scaled out anymore. On the other hand, for relevant 
examples [l^ this model dependence is negligible for the 
decay 7 r° —)> e+e “7 and we get S = 8.30 x 10“^. This can 
be rewritten in a common way as 


pNLO 


pLO 
7r°—>77 


= 0.986 X 10“% 


(60) 


Without the inclusion of the I 7 IR contribution, the 
above number would become 1.03 x 10“"^. The stated val¬ 
ues are consistent with the previous results 1.05 x 10“"^ of 
Joseph and 0.95 x 10“"^ of Mikaelin and Smith, who 
admitted that eventual numerical inaccuracy might be 
present in their result 


VII. SUMMARY 

In the preceding sections we have explored all the rel¬ 
evant NLO radiative corrections to the Dalitz decay of a 
neutral pion in the QED sector. In the direct compari¬ 
son to the earlier approach of Mikaelian and Smith Q, 
we have included into our treatment the one-photon- 
irreducible contribution. On the top of that, as an¬ 
nounced above we have enriched the vacuum polariza¬ 
tion insertion correction with the muon loop and have not 
thrown away the CJ(©) terms. The latter is connected to 
the calculation of an additional nontrivial integral. On 
the other hand, we were able to write the results in a 
more compact form even though more terms needed to 
be covered. The computational methods as well as some 
intermediate results are also provided and thus it should 
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Table I. The overall NLO correction 5{x,y) given in percent for a range of values of x and y (i.e. the Dalitz-plot corrections) 
for the process rr® —>■ e'*'e“ 7 . 


X 

y 

0.00 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

0.99 

0.01 


2.761 

2.714 

2.599 

2.449 

2.273 

2.061 

1.786 

1.402 

0.803 

-0.357 

-5.657 

0.02 


2.756 

2.720 

2.622 

2.480 

2.300 

2.073 

1.774 

1.355 

0.703 

-0.546 

-5.859 

0.03 


2.669 

2.639 

2.552 

2.419 

2.242 

2.012 

1.704 

1.267 

0.586 

-0.716 

-6.125 

0.04 


2.558 

2.531 

2.452 

2.327 

2.155 

1.925 

1.611 

1.164 

0.464 

-0.874 

-6.372 

0.05 


2.437 

2.412 

2.340 

2.221 

2.053 

1.824 

1.509 

1.054 

0.341 

-1.025 

-6.601 

0.06 


2.311 

2.288 

2.221 

2.108 

1.944 

1.717 

1.400 

0.940 

0.216 

-1.172 

-6.815 

0.07 


2.184 

2.163 

2.099 

1.990 

1.830 

1.605 

1.288 

0.824 

0.092 

-1.315 

-7.017 

0.08 


2.056 

2.036 

1.975 

1.870 

1.714 

1.491 

1.173 

0.707 

-0.033 

-1.455 

-7.211 

0.09 


1.928 

1.909 

1.851 

1.749 

1.596 

1.374 

1.057 

0.588 

-0.157 

-1.593 

-7.397 

0.10 


1.801 

1.783 

1.726 

1.628 

1.477 

1.257 

0.940 

0.469 

-0.281 

-1.729 

-7.578 

0.15 


1.170 

1.154 

1.105 

1.016 

0.874 

0.661 

0.345 

-0.131 

-0.900 

-2.394 

-8.424 

0.20 


0.546 

0.532 

0.486 

0.402 

0.266 

0.057 

-0.258 

-0.738 

-1.520 

-3.048 

-9.219 

0.25 


-0.079 

-0.092 

-0.135 

-0.217 

-0.350 

-0.556 

-0.871 

-1.355 

-2.148 

-3.704 

-9.995 

0.30 


-0.713 

-0.726 

-0.768 

-0.847 

-0.978 

-1.184 

-1.499 

-1.988 

-2.790 

-4.372 

-10.770 

0.35 


-1.366 

-1.378 

-1.419 

-1.497 

-1.627 

-1.833 

-2.149 

-2.641 

-3.454 

-5.058 

-11.558 

0.40 


-2.044 

-2.056 

-2.097 

-2.174 

-2.304 

-2.509 

-2.827 

-3.324 

-4.146 

-5.773 

-12.370 

0.45 


-2.759 

-2.771 

-2.811 

-2.887 

-3.017 

-3.222 

-3.543 

-4.044 

-4.875 

-6.525 

-13.218 

0.50 


-3.521 

-3.533 

-3.572 

-3.648 

-3.777 

-3.983 

-4.306 

-4.811 

-5.653 

-7.324 

-14.115 

0.55 


-4.344 

-4.356 

-4.395 

-4.470 

-4.599 

-4.806 

-5.130 

-5.640 

-6.492 

-8.186 

-15.076 

0.60 


-5.249 

-5.261 

-5.299 

-5.373 

-5.501 

-5.708 

-6.034 

-6.549 

-7.410 

-9.128 

-16.123 

0.65 


-6.262 

-6.273 

-6.310 

-6.383 

-6.510 

-6.717 

-7.044 

-7.563 

-8.435 

-10.177 

-17.284 

0.70 


-7.425 

-7.435 

-7.470 

-7.541 

-7.666 

-7.871 

-8.198 

-8.721 

-9.603 

-11.371 

-18.602 

0.75 


-8.802 

-8.811 

-8.844 

-8.910 

-9.031 

-9.232 

-9.558 

-10.084 

-10.976 

-12.772 

-20.143 

0.80 


-10.508 

-10.516 

-10.544 

-10.604 

-10.717 

-10.912 

-11.233 

-11.759 

-12.659 

-14.486 

-22.024 

0.85 


-12.779 

-12.784 

-12.804 

-12.851 

-12.949 

-13.129 

-13.438 

-13.958 

-14.864 

-16.724 

-24.468 

0.90 


-16.207 

-16.205 

-16.206 

-16.225 

-16.289 

-16.434 

-16.712 

-17.208 

-18.108 

-20.003 

-28.003 

0.95 


-23.167 

-23.144 

-23.084 

-23.011 

-22.960 

-22.982 

-23.140 

-23.532 

-24.360 

-26.256 

-34.451 

0.99 


-54.287 

-54.068 

-53.442 

-52.496 

-51.351 

-50.147 

-49.029 

-48.155 

-47.761 

-48.467 

-55.831 


be possible for an interested reader to trace back all the 
steps made. 

From the newly included contributions only the I 7 IR 
correction is relevant for the decay tt^ —>■ e“'"e “7 and 
should be introduced in the future analyses. Needless 
to say, the provided calculation is universal considering 
the masses of the particles involved. It can thus be also 
shown via direct calculation that if we change the masses 
of the particles in such a way that they correspond to 
the process rj —?> all the discussed corrections 

should be taken into account. In other words, both the 
muon loop as well as the 0{v'^) terms give then a non- 
negligible contribution to the overall i5(x, y). That is why 
these corrections should not be overlooked. If necessary, 
heavier charged fermions may be also introduced in the 
loops in the same way the muon loop was added. 

We believe that this work is a good starting point for 
a treatment of some other processes such as the Dalitz 
decays of 77. We have also touched on some particular 
difficulties that appear and one needs to be careful about. 
A more detailed review of this matter is beyond the scope 
of this work and will be discussed separately in the paper 
in preparation. 

Let us also say that after the complete recalculation of 
the results given in Ref. Q we have verified the formulas 


therein. The numerical accuracy of the listed values is 
also sufficient. 

The main message of the present work is the comple¬ 
tion of the list of the NLO corrections and refining of the 
expressions. All the formulas necessary for the calcula¬ 
tion of the considered correction are listed in the present 
paper in a ready-to-use form. For the eventual future 
practical use of an interested reader, we submit together 
with this text also (as ancillary files) a C++ code, which 
contains all the expressions in a well-arranged way. As a 
demonstration, the resulting program calculates the cor¬ 
rection 5(x^ y). 
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Appendix A: Bremsstrahlung matrix element squared 

For the rescaled bremsstrahlung invariant matrix element squared Tr (see ij) for the definition) we can write 

J[Tr] = J[Trc] + J[Ttu] ■ (Al) 

The explicit forms of the (from the point of view of the integration over x-y) convergent part Trc and the divergent 
part Trp of Tr are shown belovl3- The terms are already reduced (see Appendix [B] for the reduction procedure) to 
the basic integrals (see Appendix |D] for the explicit expressions) and the symmetries of the operator J (1^ and (l43l) 
were used. This means the relation Tr = Trc + T^d holds only effectively (with operator J applied). 


Trc = 16 - 4M^ 


(1 — x){3 — y) — 3xy + + 


2V‘^{X + Xy) 


(1 — a;)(l + y) — Xy\ A 


v'^M^ 


X X ) AB 


v'^\ xi B o 9 B .1 

— + 16—- - 16M^—- 


32M2 


A 


A2 




(\ — X — Xyf — AXXy — (1 — Xfy^ — 2v^{x — ?>Xy) + 


(1 + a: — Xy) E 
8iy'^{x — Xy)'^ 

(1 + a; — Xy)"^ — 2 / 2(1 — xY 


1 

AD 


8v‘^M'^{x — XyY 


(- 


{1 X — Xy)[\ + X — Xy — y(\ — x)\ \AE DE 





1 

X2 J 

AE 


\ 1 

2x J 

1 AE^ 


— + [{1 — x)y — X + Xy] —21^^} 


1 

1 1 ■ 

- 

1 

1 1 ■ 

A^E 

NPx A2_ 

A^E^ 

M4x2 A2_ 


4M4 / 1 1 

(1 + a; — Xy) \Aif CE 


(1 + 2/^)[2 — a: + a;^ — 2a;.y(l — a;)^] + 2y[l — x — xy — a:.y(l — a;^)] 


— Xy\Qx + xXy + 2(1 — a;-,,)^] + 2i^^[l + 2a; + a;.^(l + a; — a;..^)] — 

+ {y^ -y} 


AXXy{x'^ + Xy) — 2v^{x + Xy)"^ 

[(l-a;)(l + 2 /)-a;^] 


(A2) 


Tro =m4(1-x)2 ! + ?/" + 


4 1- 


1 


2x AB X \A‘^ B^ 


1 


1 


(A3) 


The integration over Xy of J[TrD] has to be done analytically. After substituting the appropriate expressions 
from (|D32I) and (ID33|) and putting the result into (I37L we find for the contribution of the divergent part to the 
bremsstrahlung correction 




2p 


O^niax 

m . 2.x 

log — + log 
A (1 “ 


- \ log(l - ^ K{x, y) \ , (A4) 


where K(x, y) is given by (ID34I) . It is apparent that the IR divergent part indeed cancels with its counterpart in the 
virtual correction 5™^{x,y). 


Note that {y —>■ —y} in Trc holds for the entire expression (in- 
eluding terms independent of y). 




































11 


Appendix B: Reduction of J terms 

In this Appendix we summarize all the necessary reductions of the J terms to the basic integrals, the results of 
which can be found in Appendix |DJ The following formulas are used to get the matrix element squared in the form 
shown in Appendix [Al 


EF 
' A 

1 

AC 

1 

'Wd 


1- ' J 

J (E + F) 

' 1 

E_ 

1 (A + C) 

(E + F) 

A' 

_E_ 

^ J 

A' 

J (A + C) 

_A_ 

1- " J 

A' 

J (B + D) 

B 


J 


M^{1 + X — Xj) 

1(1 - x)(l + y) - Xj] 
4(1 + X — Xj) 


J 


M2 [{I - x){l + y) - x^] 
8 

M2 [(1 - x)(l -y)- x^] 


L^J y^-y 


= J 


AC 


(Bl) 

(B2) 

(B3) 

(B4) 


-J v^-v 


1 


AEF 

1 

ABE 


J 


ABE^ 


(E + F) 
1 


AE 


{E-2A- 2B) 


-J 

M^x \ 2 


1 ^ 1 


r 1 1 


\ ^ 1 


- 2 J 


-2J 


[ab\ 


[ae\ 


[be\ 


AB 


1 


[E-2A- 2B) 


1 


+ J 
J 

- J 

J 


CE 


1 


M2(l + X — Xry) 


AE 


H” J 


CE 


AE 


- J 


1 


M4a;2 U 


1 


AB 


ABE 
- J 


- 2 J 


1 

AE 

1 


AE 


AE^ 

2 


M^x 


J y^-y/ 
- 2 J 

J ^ 


EE2 


AE^ 


+ {y^ -y} 


(B5) 


(B6) 


(B7) 


1 


ACE 

1 

ADE 


1 


ACEF 

1 


ADEF 


1 


(A + C) 
1 


(j 

1 

+ J 

1 

V 

[ae\ 


[ce\ 


1 


i^-A + D) \2 
4 


J 


AD 


H“ J 


AE 


M2 [(l-a;)(l+?/)-a;^] 


- J 


J 


AE 


H“ J 


CE 


DE 


M2 X — x^ — y(l — x 
2 


1 

r 1 1 


\ ^ 1 

■ 

-J 

AD 

+ J 

ae_ 

- J 

l2 


CE 


y^-y/ 


(A + C) 

1 

(E + F) 

1 


J 


1 

® (J 

1 

+ J 

1 


AEF 

M‘^{1 + X - x.y) [(1 - a:)(l + y) - a;^] V 

AE 

CE 

) 


1 


ADE 


M2(l + a; — x~f) 


+ J 
1 


1 


ADE 


BCE 
+ J 


1 


ADE 


y^-v/ 

\j 


1 M"‘(l + X — x^)[l + x — x-y — y(l — a;)] \ 2 

+ {y^ -y} 


AD 


H" J 


AE 


- J 


CE 


(B8) 


(B9) 


(BIO) 


(Bll) 
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Appendix C: Computational methods 

In this Appendix we show the approaches we used to evaluate the basic integrals listed in Appendix [P] 


1. Feynman parametrization 

With the help of the Feynman parametrization 


1 


r(ai + • • • + a„) <5 (ELi Uk - 1)<^ ' 




r(ai)---r(a„) Jq Jq + • • • + 

we can prepare, for example, the following terms for further integration 


1 2 ^ r 1 

AE ~ Jo [{E - 2A)a + 2A] 

1 pl-a 

d/3 


W2f'd„ 


A^E E(2A)(2A) ^L 


0 [2(A + ai?) + aM^x]' 

(1 — a)da 


= 2 


/o V 


da 

~2 


[{E - 2A)a + 2A]3 


1 


AA2 E^(2A) 

1 4 


= 4 


' da 


/o T 


A2A2 £;2(2A)- 


= 24 


a(l — a) da 


where we have defined r] = 21 ■ [ap + g) + aM^x + ie. Now, let us calculate the following (simplest) integral 


1 

Ie 


= 2J 


^iVo 


1 nl 

da 


f ^ 

Jo 7/2 J 47r 

dz 


dfl 


/o ^ 


1 


da 


da 


dn 


[ 2 1 ■ {ap + g) + aM^x + ie ]' 

u 21 qY 

2 da 


We have introduced 


_i (uo - |u|z + F)2 2 IIJ 0 


dloPoX dlomj (r=0) 


E + + 2uoV Jo W2a‘^ + wia + wq 


2 I {'—'JJ 2 I r/i \/i \ 1 I 2 

W2=x + = X +-x[{l-x){l-y)-Xo\ + -v Xo 


Wi = 


Wo = 


dloX dloQoX 8lQm2 (r=0) 


M2 

;2^2 


M2 


M4 


= xx^ +-x[{l-x){l + y) - x^]--v^ 


4^oTO 2 (r=0) I 2 

M4 


= —vx^. 

4 


The last integral can be evaluated as 

da 1 


_ _/ 2wo + Wi + sjwl - dwoW2 I _ J_ / p + 

Jo W2a^ + wia + wo \/w\ — dwoW 2 y2wo + wi — ^/wl^^4woW2 j \P~V^/ 

The other integrals which belong to this section can be treated in the following way 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 


(C6) 


1 

AB B{-A) 
1 (- 2 ) 


(-1) =-rda 

Jo 


1 


0 [{A + B)a — A\^ 


CE E{-2C) 

1 1 


= -2 


da 


0 [{E + 2C)a - 2Cf 


= da 


= -2 


da 


0 [2(C' + aD)+aM2(x^ - 1)]^ 


BC [{A + C)-A]B Jo [{A + B)a-B-a{A + C)y 

iHSl 


(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 
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2. Legendre polynomials and functions of the second kind 

Finally, two basic integrals can be evaluated by expanding to the Legendre functions. We can write 


B _l-p lopo - irpl cosdp 


PO 

IpI 


— cost 


a — cost 


^ ^ ^ ^oQo — 


\q\ COSUq 


On W\ |f|-cos6i, |g| b-cosOg 


(CM) 


where we have introduced a = Po/|p| and b = qo/\q\ ■ Now, consider first two Legendre polynomials and Legendre 
functions of the second kind, i.e. 


Po{x) = 1, Pi{x) = X, 


1 X 1 

Qo{x) = - log-- , Qi{x) = xQo{x) - 1. 

2 X — 1 


(C15) 


(C16) 


The numerator in (IC14p can thus be rewritten in terms of the Legendre polynomials and the denominator can be 
expanded in the following waj0 


a — cos Op _ a — I ■ p 
b — cos Oo 


b — I ■ q 

Since there is a useful integral formula for unit vectors n and Hi 


aPo(j- P^ - Pi(j- ^(2to + l)Pm(j- q^Qm{b) . (C17) 


f Aiv 

/ Pm{.H ■ni)Pm'{.H ■n 2 )dn{n) = — -— ■H 2 ), 

J (2m+1) 


(2m + 1) 

the infinite sum in (jC17ll reduces in the final result to only two terms 

I? 


(CIS) 


\P\ 


J 


47r Ini 


J ^d^ = aQoib)Po(^p- q'^ -Qiib)Pi(^p- g) = aQoib) - [p ■ g) Qi(b). (C19) 


We can differentiate the previous terms in order to get the last missing piece, since 


I - p _ \p\ a — cosOp 


B 

(Cg)2“Zo|gp (b-cosOq)'^ 


d 


a — cos 0, 


/o|g| db [|g| b — cosOq 


1 d fB 


Zo|g| db \A 


Hence 

The results can be written in the form 

J 

J 

where 


'B 

A 

B 


B 


k \m 


1 d 


J 


k|g| db 

UM'^x - 2m2) -pogo 




Po ! go 


\(M‘^x - 2m‘^) -pogo ( go 
lo\q\^ 


Qq{x) = ^ , Q[(x) = Qo(x) + xQq(x) . 


(C20) 

(C21) 

(C22) 

(C23) 

(C24) 


^ We use the hat sign to stand for the unit vector, i.e. I = I/\l \ • 
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Appendix D: Basic J terms 


In this Appendix we list the results of the basic set of 
integrals generated by acting of the operator J on the 
desired combinations of variables A, ..., F in terms of 
(1551) . First, we define a useful logarithmic function 


L(a, b) = 


1 


a>0 


\/a‘^ — b 
1 

\/o?- — b 


log 

log 


. + — b 


, — \/ a? — b 


-1 


1 ~ A / 1 ~ 


(Dl) 


and variables, in which the results have the simple forms 


= ^[{^ - x){l + y) - x^] (D2) 

V2 = ^[{l-x){l-y)-x^] (D3) 

Vq = Vl V 2 + X = -(1 + x — X^) (D4) 


p = XX J + 2xvi 
p = xxy — 2(1 — x^)vi 


K = p^ — V^XX^ 


1 


2„,2 


w = —Vq + cc + -(1 — xfy 


(D5) 

(D6) 

(D7) 

(D8) 


Co = v^{vo - 1) + P 

Cl = - 1) + xx^ 

2x 2 


C = 1- 


12 




(D9) 

(DIO) 

(Dll) 

(D12) 


Using standard integration techniques we find the fol¬ 
lowing integrals 


details see Appendix 1C 1|) we are able to calculate 
1 


AB 

1 


M'^x, 


-L{x, v'^x) 


AE 

1 




y) 


CE 

1 

1 




AE^ 

1 


M4 

4Ci 


L{uj, ly’^x-yOj) 


2C2 


J 


A^E 


M^kx M^k 

32C2 , 4Ci 


M^v'^KXry M'^k 


AE 

J 


AE 


(D18) 

(D19) 

(D20) 

(D21) 

(D22) 

(D23) 


as well as the integraljwhich did not need to be evaluated 


in the original work due to the systematic neglecting 
of the terms of order higher than 


J 


1 

16 

_A^E^_ 

“ M^k 


Ev\ 




+ ^+C 


{-ivQVi + pC)J 


AE 


(D24) 


Using the expansion to the Legendre polynomials and 
functions of the second kind (see Appendix 1C 21) , we find 


J 

J 


'B' 

= -^ + 

M'^XyC2 J 

1 


c 

2C 


A 

B ' 

A-/’ 

00 

1 

-1^ 

'1 


A\ 


A_ 

5 


where we have introduced 


Cl = 2xx^ — v'^Xj — 4:ViV2 
(2 = 2xvi — iy^(vo — x) 

C = ^v1 — v'^x-y . 


(D25) 

(D26) 


(D27) 

(D28) 

(D29) 


J 

J 


J[l] 
■ 1' 
A_ 
1 ■ 
A2 


J 

J 


'1 

A 

1 


j^L{2vi,v‘^Xy) 


16 


M^v'^x^ 


1 


2M2 

1 

■ 


L(no,x) 


(D13) 

(D14) 

(D15) 

(D16) 

(D17) 


With the help of the Feynman parametrization (for 
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We can extract the divergent parts of the integrals (ID23I) and (ID24I) through 


J 


1 


J 


A^E 

1 


1 


M'^x 

1 

4 


J 


1 

1 


16^0 , 4a 


M^v'^k M'^k 


AE 


A2 


16 


1 + 


1 




4ui + + 


3xj^q 

2k 


M^k 


(4voVi 


AE 


(D30) 

(D31) 


The above formulas have a very convenient form and are to be substituted into ___ 

The divergent integrals alone have to be integrated over Xj analytically. The calculation is done in detail in Ref. @ 
and the results can be written in a simple form 


J 


J 


1 

1 

AB 


da;.v = 


16 


^ 1/2^4 


22,max 


1 m 

log—+ log 

A (l-x)(l+y)J 


da;.^ = — ——;—^ 

M‘^x/3 


O^max 

m . 2.x 

log — + log 
A (1 ~ 


log( 7 ) 




K{x,y ), 


where 


K(x,y) = 


2 log 


y+ /3 

2/3 


log — 

X 


log(7) - Li2 


i{y- P) 


. y + P _ 


Li2 


y-13 
.7(y+ /3). 


(D32) 

(D33) 

(D34) 


These terms are to be used to evaluate f J[TrD] da;.y . 
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